Statistical modeling of the eye for multimodal treatment planning for external beam radiation therapy of intraocular tumors

Ocular anatomy and radiation-associated toxicities provide unique challenges for external beam radiation therapy. For treatment planning, precise modeling of organs at risk and tumor volume are crucial. Development of a precise eye model and automatic adaptation of this model to patients' anatomy remain problematic because of organ shape variability. This work introduces the application of a 3-dimensional (3D) statistical shape model as a novel method for precise eye modeling for external beam radiation therapy of intraocular tumors.

Enhancing the Visualization of the Peripheral Retina with Wide Field-of-View Optical Coherence Tomography

The goal of my Ph.D. thesis is to enhance the visualization of the peripheral retina using wide-field optical coherence tomography (OCT) in a clinical setting.

OCT has gain widespread adoption in clinical ophthalmology due to its ability to visualize the diseases of the macula and central retina in three-dimensions, however, clinical OCT has a limited field-of-view of 300. There has been increasing interest to obtain high-resolution images outside of this narrow field-of-view, because three-dimensional imaging of the peripheral retina may prove to be important in the early detection of neurodegenerative diseases, such as Alzheimer's and dementia, and the monitoring of known ocular diseases, such as diabetic retinopathy, retinal vein occlusions, and choroid masses.

Before attempting to build a wide-field OCT system, we need to better understand the peripheral optics of the human eye. Shack-Hartmann wavefront sensors are commonly used tools for measuring the optical imperfections of the eye, but their acquisition speed is limited by their underlying camera hardware. The first aim of my thesis research is to create a fast method of ocular wavefront sensing such that we can measure the wavefront aberrations at numerous points across a wide visual field. In order to address aim one, we will develop a sparse Zernike reconstruction technique (SPARZER) that will enable Shack-Hartmann wavefront sensors to use as little as 1/10th of the data that would normally be required for an accurate wavefront reading. If less data needs to be acquired, then we can increase the speed at which wavefronts can be recorded.

For my second aim, we will create a sophisticated optical model that reproduces the measured aberrations of the human eye. If we know how the average eye's optics distort light, then we can engineer ophthalmic imaging systems that preemptively cancel inherent ocular aberrations. This invention will help the retinal imaging community to design systems that are capable of acquiring high resolution images across a wide visual field. The proposed model eye is also of interest to the field of vision science as it aids in the study of how anatomy affects visual performance in the peripheral retina.

Using the optical model from aim two, we will design and reduce to practice a clinical OCT system that is capable of imaging a large (800) field-of-view with enhanced visualization of the peripheral retina. A key aspect of this third and final aim is to make the imaging system compatible with standard clinical practices. To this end, we will incorporate sensorless adaptive optics in order to correct the inter- and intra- patient variability in ophthalmic aberrations. Sensorless adaptive optics will improve both the brightness (signal) and clarity (resolution) of features in the peripheral retina without affecting the size of the imaging system.

The proposed work should not only be a noteworthy contribution to the ophthalmic and engineering communities, but it should strengthen our existing collaborations with the Duke Eye Center by advancing their capability to diagnose pathologies of the peripheral retinal.

Modeling the Possible Influences of Eye Movements on the Refinement of Cortical Direction Selectivity

The second-order statistics of neural activity was examined in a model of the cat LGN and V1 during free-viewing of natural images. In the model, the specific patterns of thalamocortical activity required for a Bebbian maturation of direction-selective cells in VI were found during the periods of visual fixation, when small eye movements occurred, but not when natural images were examined in the absence of fixational eye movements. In addition, simulations of stroboscopic reming that replicated the abnormal pattern of eye movements observed in kittens chronically exposed to stroboscopic illumination produced results consistent with the reported loss of direction selectivity and preservation of orientation selectivity. These results suggest the involvement of the oculomotor activity of visual fixation in the maturation of cortical direction selectivity.

Optical modeling of aspherical lenses: a human eye study

We have been developing a computational code to project optical lenses, with low aberration effects. Our main interest is model the human eye, particularly, project special corrective lenses. As the lens shape is the focus of the optimization, we have coupled a ray tracing method with Monte Carlo techniques. The initial results indicated that the algorithm must be improved in terms of resolution and reliability.

Modeling the role of fixational eye movements in real-world scenes

Our eyes never remain still. Even when we stare at a fixed point, small involuntary movements take place in our eyes in an imperceptible manner. Researchers agree on the presence of three main contributions to eye movements when we fix the gaze: microsaccades, drifts and tremor. These small movements carry the image across the retina stimulating the photoreceptors and thus avoiding fading. Nowadays it is commonly accepted that these movements can improve the discrimination performance of the retina. In this paper, several retina models with and without fixational eye movements were implemented by mean of RetinaStudio tool to test the feasibility of these models to be incorporated in future neuroprostheses. For this purpose each retina model has been stimulated with natural scene images in two experiments. Results are discussed from the point of view of a neuroprosthesis development.

Three-dimensional modeling of the human eye based on magnetic resonance imaging

PURPOSE. A methodology for noninvasively characterizing the three-dimensional (3-D) shape of the complete human eye is not currently available for research into ocular diseases that have a structural substrate, such as myopia. A novel application of a magnetic resonance imaging (MRI) acquisition and analysis technique is presented that, for the first time, allows the 3-D shape of the eye to be investigated fully. METHODS. The technique involves the acquisition of a T2-weighted MRI, which is optimized to reveal the fluid-filled chambers of the eye. Automatic segmentation and meshing algorithms generate a 3-D surface model, which can be shaded with morphologic parameters such as distance from the posterior corneal pole and deviation from sphericity. Full details of the method are illustrated with data from 14 eyes of seven individuals. The spatial accuracy of the calculated models is demonstrated by comparing the MRI-derived axial lengths with values measured in the same eyes using interferometry. RESULTS. The color-coded eye models showed substantial variation in the absolute size of the 14 eyes. Variations in the sphericity of the eyes were also evident, with some appearing approximately spherical whereas others were clearly oblate and one was slightly prolate. Nasal-temporal asymmetries were noted in some subjects. CONCLUSIONS. The MRI acquisition and analysis technique allows a novel way of examining 3-D ocular shape. The ability to stratify and analyze eye shape, ocular volume, and sphericity will further extend the understanding of which specific biometric parameters predispose emmetropic children subsequently to develop myopia. Copyright © Association for Research in Vision and Ophthalmology.

Structural Integrity of the Greek Key Motif in beta gamma-Crystallins Is Vital for Central Eye Lens Transparency

Background: We highlight an unrecognized physiological role for the Greek key motif, an evolutionarily conserved super-secondary structural topology of the beta gamma-crystallins. These proteins constitute the bulk of the human eye lens, packed at very high concentrations in a compact, globular, short-range order, generating transparency. Congenital cataract (affecting 400,000 newborns yearly worldwide), associated with 54 mutations in beta gamma-crystallins, occurs in two major phenotypes nuclear cataract, which blocks the central visual axis, hampering the development of the growing eye and demanding earliest intervention, and the milder peripheral progressive cataract where surgery can wait. In order to understand this phenotypic dichotomy at the molecular level, we have studied the structural and aggregation features of representative mutations. Methods: Wild type and several representative mutant proteins were cloned, expressed and purified and their secondary and tertiary structural details, as well as structural stability, were compared in solution, using spectroscopy. Their tendencies to aggregate in vitro and in cellulo were also compared. In addition, we analyzed their structural differences by molecular modeling in silico. Results: Based on their properties, mutants are seen to fall into two classes. Mutants A36P, L45PL54P, R140X, and G165fs display lowered solubility and structural stability, expose several buried residues to the surface, aggregate in vitro and in cellulo, and disturb/distort the Greek key motif. And they are associated with nuclear cataract. In contrast, mutants P24T and R77S, associated with peripheral cataract, behave quite similar to the wild type molecule, and do not affect the Greek key topology. Conclusion: When a mutation distorts even one of the four Greek key motifs, the protein readily self-aggregates and precipitates, consistent with the phenotype of nuclear cataract, while mutations not affecting the motif display `native state aggregation', leading to peripheral cataract, thus offering a protein structural rationale for the cataract phenotypic dichotomy ``distort motif, lose central vision''.

A common stochastic accumulator with effector-dependent noise can explain eye-hand coordination

The computational architecture that enables the flexible coupling between otherwise independent eye and hand effector systems is not understood. By using a drift diffusion framework, in which variability of the reaction time (RT) distribution scales with mean RT, we tested the ability of a common stochastic accumulator to explain eye-hand coordination. Using a combination of behavior, computational modeling and electromyography, we show how a single stochastic accumulator to threshold, followed by noisy effector-dependent delays, explains eye-hand RT distributions and their correlation, while an alternate independent, interactive eye and hand accumulator model does not. Interestingly, the common accumulator model did not explain the RT distributions of the same subjects when they made eye and hand movements in isolation. Taken together, these data suggest that a dedicated circuit underlies coordinated eye-hand planning.

Analog VLSI-based, neuromorphic sensorimotor systems: modeling the primate oculomotor system

Using neuromorphic analog VLSI techniques for modeling large neural systems has several advantages over software techniques. By designing massively-parallel analog circuit arrays which are ubiquitous in neural systems, analog VLSI models are extremely fast, particularly when local interactions are important in the computation. While analog VLSI circuits are not as flexible as software methods, the constraints posed by this approach are often very similar to the constraints faced by biological systems. As a result, these constraints can offer many insights into the solutions found by evolution. This dissertation describes a hardware modeling effort to mimic the primate oculomotor system which requires both fast sensory processing and fast motor control. A one-dimensional hardware model of the primate eye has been built which simulates the physical dynamics of the biological system. It is driven by analog VLSI circuits mimicking brainstem and cortical circuits that control eye movements. In this framework, a visually-triggered saccadic system is demonstrated which generates averaging saccades. In addition, an auditory localization system, based on the neural circuits of the barn owl, is used to trigger saccades to acoustic targets in parallel with visual targets. Two different types of learning are also demonstrated on the saccadic system using floating-gate technology allowing the non-volatile storage of analog parameters directly on the chip. Finally, a model of visual attention is used to select and track moving targets against textured backgrounds, driving both saccadic and smooth pursuit eye movements to maintain the image of the target in the center of the field of view. This system represents one of the few efforts in this field to integrate both neuromorphic sensory processing and motor control in a closed-loop fashion.

Computational modeling and psychophysics in low- and mid-level vision

This thesis addresses a series of topics related to the question of how people find the foreground objects from complex scenes. With both computer vision modeling, as well as psychophysical analyses, we explore the computational principles for low- and mid-level vision.

We first explore the computational methods of generating saliency maps from images and image sequences. We propose an extremely fast algorithm called Image Signature that detects the locations in the image that attract human eye gazes. With a series of experimental validations based on human behavioral data collected from various psychophysical experiments, we conclude that the Image Signature and its spatial-temporal extension, the Phase Discrepancy, are among the most accurate algorithms for saliency detection under various conditions.

In the second part, we bridge the gap between fixation prediction and salient object segmentation with two efforts. First, we propose a new dataset that contains both fixation and object segmentation information. By simultaneously presenting the two types of human data in the same dataset, we are able to analyze their intrinsic connection, as well as understanding the drawbacks of today’s “standard” but inappropriately labeled salient object segmentation dataset. Second, we also propose an algorithm of salient object segmentation. Based on our novel discoveries on the connections of fixation data and salient object segmentation data, our model significantly outperforms all existing models on all 3 datasets with large margins.

In the third part of the thesis, we discuss topics around the human factors of boundary analysis. Closely related to salient object segmentation, boundary analysis focuses on delimiting the local contours of an object. We identify the potential pitfalls of algorithm evaluation for the problem of boundary detection. Our analysis indicates that today’s popular boundary detection datasets contain significant level of noise, which may severely influence the benchmarking results. To give further insights on the labeling process, we propose a model to characterize the principles of the human factors during the labeling process.

The analyses reported in this thesis offer new perspectives to a series of interrelating issues in low- and mid-level vision. It gives warning signs to some of today’s “standard” procedures, while proposing new directions to encourage future research.

Cooperation of Mtmr8 with PI3K Regulates Actin Filament Modeling and Muscle Development in Zebrafish

Background: It has been shown that mutations in at least four myotubularin family genes (MTM1, MTMR1, 2 and 13) are causative for human neuromuscular disorders. However, the pathway and regulative mechanism remain unknown. Methodology/Principal Findings: Here, we reported a new role for Mtmr8 in neuromuscular development of zebrafish. Firstly, we cloned and characterized zebrafish Mtmr8, and revealed the expression pattern predominantly in the eye field and somites during early somitogenesis. Using morpholino knockdown, then, we observed that loss-of-function of Mtmr8 led to defects in somitogenesis. Subsequently, the possible underlying mechanism and signal pathway were examined. We first checked the Akt phosphorylation, and observed an increase of Akt phosphorylation in the morphant embryos. Furthermore, we studied the PH/G domain function within Mtmr8. Although the PH/G domain deletion by itself did not result in embryonic defect, addition of PI3K inhibitor LY294002 did give a defective phenotype in the PH/G deletion morphants, indicating that the PH/G domain was essential for Mtmr8's function. Moreover, we investigated the cooperation of Mtmr8 with PI3K in actin filament modeling and muscle development, and found that both Mtmr8-MO1 and Mtmr8-MO2+LY294002 led to the disorganization of the actin cytoskeleton. In addition, we revealed a possible participation of Mtmr8 in the Hedgehog pathway, and cell transplantation experiments showed that Mtmr8 worked in a non-cell autonomous manner in actin modeling. Conclusion/Significance: The above data indicate that a conserved functional cooperation of Mtmr8 with PI3K regulates actin filament modeling and muscle development in zebrafish, and reveal a possible participation of Mtmr8 in the Hedgehog pathway. Therefore, this work provides a new clue to study the physiological function of MTM family members.

Heritability of peripheral refraction in Chinese children and adolescents: the Guangzhou Twin Eye study.

PURPOSE:

To estimate the heritability of peripheral refraction in Chinese children and adolescents.

METHODS:

The authors examined 72 monozygotic (MZ) twins and 48 dizygotic (DZ) twins aged 8 to 20 years from a population-based twin registry. Temporal and nasal peripheral refraction, each 40° from the visual axis, and axial refraction were measured using an autorefractor. Relative peripheral refractive error (RPRE) was defined as the peripheral refraction minus the axial refraction. Heritability was assessed by structural equation modeling after adjustment for age and sex.

RESULTS:

The mean and SD of temporal refraction (T(40)), nasal refraction (N(40)), RPRE-T(40), RPRE-N(40), and T(40)-N(40) asymmetry were -0.27 ± 2.0 D, 0.36 ± 2.19 D, 1.18 ± 1.39 D, 1.80 ± 1.69 D, and -0.62 ± 1.58 D, respectively. The intraclass correlations for T(40) refraction, N(40) refraction, RPRE-T(40), RPRE-N(40), and T(40)-N(40) asymmetry were 0.87, 0.83, 0.65, 0.74, and 0.58 for MZ pairs and 0.49, 0.42, 0.30, 0.41, and 0.32 for DZ pairs, respectively. A model with additive genetic and unique environmental effects was the most parsimonious, with heritability values estimated as 0.84, 0.76, 0.63, 0.70, and 0.55, respectively, for the peripheral refractive parameters.

CONCLUSIONS:

Additive genetic effects appear to explain most of the variance in peripheral refraction and relative peripheral refraction when adjusting for the effects of axial refraction.

Contribution of genetic and environmental effects on lens thickness: the Guangzhou Twin Eye study.

PURPOSE:

This study investigated the heritability of lens thickness (LT) and relative lens thickness (LT/axial length, rLT) measured by Lenstar among Chinese children and adolescents in the Guangzhou Twin Eye study.

METHODS:

Twins aged 8 to 22 years were enrolled from the Guangzhou Twin Registry. A series of LT and axial length (AL) measurements using the Lenstar were taken for each twin. Zygosity was confirmed by genotyping in all same-sex twin pairs. Heritability was assessed by structural variance component genetic modeling, after adjustment for age and sex with the Mx program.

RESULTS:

Seven hundred sixty-eight twin pairs (482 monozygotic [MZ] and 286 dizygotic [DZ] twins) were available for data analysis. The mean (standard deviation) LT and rLT were 3.45 (0.18) mm and 0.142 (0.01), respectively. The intraclass correlation coefficients (ICCs) for LT were 0.90 for the MZ and 0.39 for the DZ twins; and those for rLT were 0.90 for the MZ and 0.40 for the DZ twins, respectively. The best-fitting model yielded 89.5% (95% CI: 87.8%-91.0%) of additive genetic effects and 10.5% (95% CI: 9.0%-12.2%) of unique environmental effects for LT, and 89.3% (95% CI: 89.2%-89.3%) of additive genetic effects and 10.7% (95% CI: 10.7%-11.4%) of unique environmental effects for rLT.

CONCLUSIONS:

This study confirms that the LT in young healthy subjects may be mainly affected by additive genetic factors. High heritability remains even when the data are corrected for the influence of AL with the use of rLT.

Metagenomic systems biology: frameworks for modeling and characterizing the gut microbiome

Thesis (Ph.D.)--University of Washington, 2014